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Abstract 

The  research  described  in  this  report  represents  an  important  step  forward  in  the 
development  of  seabed  classification  technology.  First,  it  represents  the 
application  of  the  QTC  VIEW  technology  to  a  new  seabed  environment.  The  data 
sets  analyzed  here  involve  reef  seabeds  which  have  not  previously  been  studied. 
Second,  they  represent  a  consistent  effort  to  compare  the  results  obtained  by 
using  two  different  frequencies  of  sounder  over  the  same  seabed  type.  The 
results  suggest  that  both  200  kHz  and  24  kHz  sounders  can  be  used  to  classify  the 
four  bottom  types  involved  in  the  study.  However,  it  is  clear  that  multispectral 
seabed  classification  holds  significant  potential  in  situations  where,  for  example,  a 
gravel  bedform  is  overlain  by  a  thin  veneer  of  mud. 


Resume 

Les  recherches  que  nous  ddcrivons  dans  ce  rapport  representent  une  etape  en 
avant  significative  dans  le  developpement  de  la  technologic  pour  la  classification 
des  lits  de  la  mer.  D'abord,  elles  representent  I'application  de  la  technologic 
QTC  VIEW  a  un  nouvel  environnement  de  lits  de  la  mer.  Les  ensembles  de 
donnees  analyses  ici  comprennent  des  lits  de  la  mer  a  recifs,  que  Ton  n'a  pas 
etudies  auparavant.  Ensuite,  nos  recherches  representent  un  effort  constant  de 
comparer  les  resultats  obtenus  en  utilisant  deux  frequences  de  sondeur  differentes 
au-dessus  d'un  meme  type  de  lit  de  la  mer.  Les  resultats  suggerent  que  les 
sondeurs  a  200  kHz  et  a  24  kHz  peuvent,  tous  les  deux,  servir  a  classifier  les 
quatre  types  de  fond  impliques  dans  I'etude.  Neanmoins,  c'est  claire  que  la 
classification  multispectrale  de  lits  de  la  mer  montre  du  potentiel  significatif  dans 
les  situations  ou,  par  exemple,  un  lit  de  gravier  est  reconvert  d'une  couche  mince 
de  boue. 


QUE5TER 

JflnCEnT 


March  1996 


Page  V 


Florida  Keys  Seabed  Classification  Study 


SC  96-002 


1 .  Introduction 

This  report  is  supplied  under  the  terms  of  a  contract  managed  by  the  Esquimalt 
Defence  Research  Detachment  (EDRD)  of  the  Defence  Research  Establishment 
Atlantic  (DREA),  and  supplied  by  Supply  and  Services  Canada.  The  contract, 
number  W7708-5-0208/01-XSA,  involving  research  into  depth  dependence  of 
Quester  Tangent  Corporation's  (QTC)  seabed  classification  system  and  seabed 
classification  from  echo  soundings  acquired  in  the  Dry  Tortugas  area  near  Key 
West,  Florida,  was  completed  by  March  31,  1996.  This  report  describes  research 
background,  methodology,  and  results  of  Part  2,  Florida  Keys  data. 

The  research  described  in  this  report  represents  an  important  step  forward  in  the 
development  of  seabed  classification  technology.  First,  it  represents  the 
application  of  the  QTC  VIEW  technology  to  a  new  seabed  environment.  The  data 
sets  analyzed  here  involve  reef  seabeds  which  have  not  previously  been  studied. 
Second,  they  represent  a  consistent  effort  to  compare  the  results  obtained  by 
using  two  different  frequencies  of  sounder  over  the  same  seabed  type.  The 
results  suggest  that  both  200  kHz  and  24  kHz  sounders  can  be  used  to  classify  the 
four  bottom  types  involved  in  the  study.  However,  it  is  clear  that  multispectral 
seabed  classification  holds  significant  potential  in  situations  where,  for  example,  a 
gravel  bedform  is  overlain  by  a  thin  veneer  of  mud. 

Section  2  of  this  report  describes  the  geological  setting  and  ground  truth  from  the 
Florida  Keys,  Dry  Tortugas  study  area  (see  Figure  1.1).  A  3'  x  3'  grid  was 
surveyed  at  1/2'  spacing  resulting  in  the  acquisition  of  over  70  line  km  of  acoustic 
data.  A  wide  variety  of  bottom  information  was  collected  for  the  work  using 
several  vessels,  and  remote  and  in-situ  measurement  techniques.  While  the 
ground  truth  is  not  evenly  spaced  over  the  study  area,  the  ground  truth 
information  is  among  the  most  detailed  available  for  this  type  of  work  to  date. 

Section  3  describes  the  seabed  classification  techniques  used  to  process  the 
acoustic  data,  and  Section  4  describes  classification  results.  Sufficient  detail  of 
results  is  presented  on  a  line-by-line  basis  to  allow  the  reader  to  draw  conclusions 
as  to  the  classification  accuracy  obtained  in  the  study. 

Finally,  Section  5  presents  conclusions  and  suggestions  for  further  research. 
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2 .  Geological  Interpretation  and  Ground  Truthing 

This  section  describes  the  geological  setting  of  the  study  area  and  ground  truthing 
available  to  confirm  the  seabed  classification  produced  by  the  QTC  seabed 
classification  system. 


2.1  Geological  Setting 

The  study  area  (  24°35'N,  82“50'W)  is  located  just  south  of  the  Dry  Tortugas 
area  at  the  western  leading  edge  of  the  Florida  Keys  (Figure  1.1).  The  Keys  are 
a  series  of  coral  sand  cays  and  exposed  limestone  reef  flats  extending  westward 
from  the  Atlantic  coast  of  Florida  along  the  southern  edge  of  Rorida  Bay.  The 
cays  of  the  Dry  Tortugas  are  sub-aerially  exposed  lithified  sand  shoals.  The  cays 
range  from  2  to  10  metres  above  sea  level.  Much  of  the  coastline  consists  of 
mangrove  swamps. 

The  study  area  is  located  south  of  Garden  Key  in  water  depths  ranging  from  5m 
to  30m.  The  area  is  underlain  by  limestone  bedrock  which  has  been  dated  at 
125,000  BP.  At  6000  years  BP  the  sea  level  was  relatively  lower  and  a  series  of 
fringing  reefs  were  formed  in  the  southwest  of  the  survey  area.  A  relatively 
quick  rise  in  sea  level  coupled  with  a  high  influx  of  sediment  likely  proved  hostile 
to  reef  growth  and  may  explain  why  the  reefs  failed  to  thrive.  Remnants  of  these 
fringing  reefs  can  be  seen  in  the  survey  area.  The  largest  of  the  three  reef  flat 
features  trends  northeast-southwest  and  is  approximately  1  km  wide.  The 
northern  and  southern  extent  of  the  reef  flats  are  beyond  the  limits  of  the  survey 
area.  The  surficial  sediments  are  biogenic  and  are  generally  silts  and  clayey-silts 
with  minor  sand  and  gravel  deposits.  Sediment  distribution  is  likely  controlled 
by  storm  deposition  and  modified  by  tidal  induced  circulation. 


2.2  Ground  Truthing  Data  and  Processing 
2.2.1  Ground  Truthing 


Data  from  58  sample  sites  were  used  as  the  basis  for  the  description  of  the 
surficial  geology  (Table  2.1).  Full  descriptions  of  the  sample  site  data  can  be 
found  in  Appendix  A.  The  sample  sites  were  positioned  using  GPS.  Sample 
distribution  was  patchy  with  the  majority  of  samples  collected  within  a  1  km 
radius  in  the  northwest  comer  of  the  survey  area  (Figure  2.1).  Data  were 
collected  using  three  vessels,  the  Planet,  the  Pelican  and  the  Seward  Johnson. 

Sampling  methods  included  gravity  and  diver  core  sampling,  and  cone 
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penetrometer  and  diver-held  shear  vane  measurements.  In-situ  pressure  and 
shear  wave  velocities  were  measured  at  selected  sites.  Sidescan  sonar  and  3.5  kHz 
sub-bottom  profiling  surveys  were  completed  in  addition  to  a  towed- video  survey 
and  still  camera  photographs.  There  was  no  towed-body  positioning  for  the 
video  survey. 

Selected  samples  were  analyzed  to  provide  mean  grainsize  and  percent  gravel, 
sand,  silt  and  clay.  Physical  property  analysis  was  used  to  provide  wet  bulk 
density,  grain  density,  water  content,  void  ratio  and  porosity.  In  some  cases  the 
physical  property  analysis  along  with  the  velocity  measurements  provided  enough 
information  to  calculate  the  acoustic  reflection  coefficient. 


2.2.2  Map  Preparation 

A  base  map  for  the  area  was  prepared  in  a  GIS  (Mapinfo).  The  reef  areas  plotted 
on  the  base  map  are  from  interpreted-sidescan  imagery  made  available  to  QTC. 
The  interpretation  was  provided  on  a  facsimile  copy;  the  reef  locations  were 
digitized  from  the  facsimile  and  imported  into  the  GIS.  The  reef  locations  must 
therefore  be  considered  approximate.  The  sample  locations  and  information 
were  imported  into  the  GIS  from  the  Excel  spreadsheets  which  were  provided. 

The  bathymetry  was  produced  from  the  acoustic  and  positioning  information 
provided  by  the  QTC  ISAH-S  system  (see  Sections  3.2  -  Data  Acquisition  and  3.3 
-  Data  Pre-Processing).  The  boat  track  is  plotted  in  Figure  2.2.  The  GPS 
positioning  (unfortunately  generally  not  differentially  corrected)  was  processed 
using  QTCs  HYPS  processing  software.  Water  depths  were  calculated  using 
QTCs  Basic  Waveform  Processing  software.  The  raw  water  depths,  uncorrected 
for  tide  and  draft,  and  not  tied  to  any  datum  were  imported  into  surface  modeling 
software  (Surfer)  and  gridded  using  triangulation  with  linear  interpolation. 
Contours  were  produced  and  smoothed  using  a  spline  function.  These  contours 
were  then  imported  into  the  GIS.  Unfortunately  the  positioning  is  not  precise 
enough  to  produce  an  accurate  bathymetry  (QTC  software  can  processes  the 
ISAH-S  acoustic  data  to  a  depth  accuracy  of  better  than  10  cm  assuming  the  sound 
speed  profile  is  known);  however,  the  relative  bathymetry  is  adequate  for  the 
purposes  of  this  paper. 
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Table  2. 1 :  Compilation  of  the  sediment  analysis  and  in-situ  measurements.  The 
grainsize  results  are  from  the  uppermost  section  of  the  cores. 


Sample 


I? 


EiCTiaTOiai'a 


147IKW-PE-GC 


167  KW-PE-GC 


178  KW-PE-GC 


203  KW-PE-GC 


210  KW-PE-GC 


213  KW-PE-GC 


218  KW-PE-GC 


222  KW-PE-GC 


224  KW-PE-GC 


225  KW-PE-GC 


227  KW-PE-GC 


208  KW-PE-GC 


220  KW-PE-GC 


179  KW-PE-GC 


285  KW-PE-GC 


113  KW-PL-DC 


173  KW-PL-DC 


192  KW-PL-DC 


198  KW-PL-DC 


215  KW-PL-DC 


223  KW-PL-DC 


244  KW-PL-DC 


263  KW-PL-DC 


236  KW-PL-DC-IA 


241  KW-PL-DC-IA 


242  KW-PL-DC-IA 


251  KW-PL-DC-IA 


253  KW-PL-DC-IA 


KW-PL-DC-IA 


257  KW-PL-DC-IA 


265  KW-PL-DC-IA 


266  KW-PL-DC-IA 


2741  KW-PL-DC-IA 


KW-PUX:-IA 


284  KW-PL-DC-IA 


286  KW-PL-DC-IA 


288  KW-PL-DC-IA 


290  KW-PL-DC-IA 


248  KW-SJ-GS-XB 


KW-SJ-GC-XB 


254  KW-SJ-GS-XB 


255 


244  KW-SJ-PCT 


265  KW-SJ-PCT 


267  KW-SJ-PCT 


297  KW-SJ-PCT 


319  KW-SJ-PCT 


Refl. 


Coeff. 


61.26  0.277 


311 

KW-SJ-GC 

326 

KW-SJ-GC 

m 
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2.3  Results 

2.3.1  Relative  Bathymetry 

The  relative  bathymetry  over  the  survey  area  is  shown  in  Figure  2.3.  The  seabed 
exhibits  approximately  11m  of  relief  over  the  survey  area.  There  is  a  general 
shallowing  trend  from  southeast  to  northwest;  the  seabed  is  relatively  flat  in  the 
northern  section  and  between  the  reef  areas.  The  most  prominent  bathymetric 
features  are  pinnacles  which  occur  in  the  southern  section.  A  4m  pinnacle  is 
associated  with  the  southern  end  of  a  reef  flat  area.  This  feature  is  about  100m  x 
400m  in  size  with  the  long  axis  oriented  parallel  to  the  central  axis  of  the  flat  reef 
area.  A  series  of  three  smaller  pinnacles  of  approximately  3m  relief  are  aligned 
parallel  to  the  central  axis  of  the  main  reef  and  appear  to  be  associated  with  the 
small  reef  flat  area  to  the  southeast.  Unfortunately  this  could  not  be  confirmed  as 
the  pinnacles  occur  beyond  the  southern  limit  of  the  sidescan  data  which  were 
provided.  It  is  tempting  to  interpret  these  pinnacles  as  modem  reef  colonies 
because  old  reef  flats  often  provide  a  stable  platform  conducive  to  reef  growth. 


2.3.2  Sediment  Texture 


Figure  2.4  shows  sediment  type  for  selected  samples.  Two  samples  at  the 
southern  end  of  the  central  reef  flat  area  were  described;  one  was  silt  and  the 
other  was  coral,  rock.  The  coral,  rock  sample  is  consistent  with  the  interpreted 
reef  flat  location.  The  sample  of  silt  may  have  been  collected  from  deposits  of 
fine  sediment  which  accumulate  in  the  interstitial  voids  in  a  reef  flat. 
Alternatively,  the  coral  to  silt  transition  between  the  two  samples  could  indicate 
that  the  reef  flat  boundary  as  interpreted  from  the  sidescan  sonar  is  inaccurate. 

From  sample  248,  located  at  the  north  end  of  the  central  reef  flat,  the  sediments 
become  progressively  finer  towards  the  northwest,  ranging  from  medium  sand 
through  silty  sand,  sandy  silt  and  clayey  silt,  the  more  sandy  samples  were 
collected  from  the  broad  shallow  area  located  between  the  north  and  central  reef 
flats.  This  area  may  represent  a  higher  energy  environment. 


Most  samples  show  little  gravel  in  the  sediment.  Apart  from  one  site 
corresponding  to  a  bathymetric  high  at  the  north  end  of  the  central  reef  flat  the 
gravel  content  in  the  sediments  seems  to  be  highest  at  the  head  of  the  northeast 
-trending  trough  at  the  west  of  the  survey  area.  The  percentage  of  sand  in  the 
seabed  sediments  is  at  a  maximum  of  90%  to  the  north  of  the  central  reef  flat. 

The  average  percentage  of  sand  in  the  rest  of  the  samples  was  about  24%.  There 
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is  a  trend  towards  a  decrease  in  the  percentage  of  sand  towards  the  shallower 
water  depths  to  the  northwest.  The  percentage  of  silt  increases  to  the  northwest. 
The  area  adjacent  to  the  western  edge  of  the  northernmost  reef  flat  contains  the 
highest  percentage  of  clay.  Figures  2.5  -  2.8  show  the  percentages  of  gravel, 
sand,  silt  and  clay. 

Table  2.2  is  a  compilation  of  the  percentage  of  gravel  size  fraction  in  the 
sediments  to  60  cm  below  the  seabed.  This  table  shows  a  broad  increase  in  the 
percentage  of  gravel  with  depth.  Note  the  significant  increase  in  the  amount  of 
gravel  at  40  cm  to  60  cm  below  the  seabed. 


Sample 

2  cm 

10  cm 

20  cm 

30  cm 

60  cm 

2-30  cm 

147 

0.00 

■ElEi 

10.00 

2.20 

3.80 

4.50 

3.08 

167 

0.50 

3.20 

2.30 

4.50 

3.90 

7.90 

4.70 

2.63 

178 

2.30 

0.67 

0.68 

IHSI 

1.31 

3.06 

1.96 

1.29 

203 

4.40 

25.30 

13.80 

10.00 

8.40 

5.10 

0.08 

13.38 

210 

0.58 

1.79 

1.55 

1.25 

6.46 

3.19 

4.37 

1.29 

224 

0.18 

0.50 

0.81 

1.59 

1.18 

1.39 

2.83 

0.77 

113 

0.03 

0.06 

0.36 

0.54 

0.25 

173 

0.05 

0.39 

2.81 

0.76 

1.00 

192 

0.69 

0.00 

0.00 

0.00 

0.17 

198 

0.00 

0.10 

0.16 

0.70 

0.24 

244 

0.00 

0.18 

0.00 

0.38 

0.14 

Table  2.2:  Percent  gravel  at  increasing  core  depth  for  selected  sites. 

2.3.3  Video  Imagery 

Approximately  45  minutes  of  video  imagery  was  analyzed.  One  video  transect 
across  the  north  area  of  the  central  reef  flat  showed  a  relatively  flat  seabed  with 
little  relief.  The  surficial  sediments  consist  of  a  veneer  of  fine  to  medium  sand 
over  a  reef  flat.  Small  cobbles  of  apparent  reef  debris  lay  on  the  surface  of  the 
reef  flat  in  some  areas.  The  reef  flat  exhibits  a  rough  surface  texture  with 
cavities  on  the  order  of  10  cm  in  diameter. 

A  second  video  transect  from  an  area  to  the  west  of  the  northern  reef  flat  was 
also  analyzed.  The  seabed  is  again  relatively  flat  with  little  surface  roughness. 
The  surficial  sediments  appear  to  consist  of  very  fine  grained  material.  A 
preponderance  of  flora  on  the  seabed  is  evident. 
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2.4  Discussion 

The  very  patchy  nature  of  the  sampling  makes  a  detailed  analysis  of  the  sediment 
distribution  difficult.  The  dominant  component  of  the  geology  within  the  survey 
area  are  the  reef  flats  which  have  likely  greatly  influenced  the  Holocene  sediment 
deposition.  The  reef  flats  are  essentially  sediment  free  except  for  some  reef 
debris  which  cannot  be  transported  except  under  extreme  conditions. 

In  a  fringing  reef  environment  material  is  eroded  from  the  seaward  edge  of  the 
reef  and  transported  over  the  reef  crest  and  deposited  in  backreef  or  lagoonal 
areas.  Tide-induced  circulation  patterns  in  backreef  areas  tend  to  transport 
material  in  a  shore-parallel  direction  until  the  material  is  flushed  from  the  system 
through  a  reef  passage  or  by  deposition  on  beaches  through  constructional  wave 
activity.  In  pre-Holocene  times  when  sea  level  was  lower  a  higher  energy 
environment  existed  and  may  account  for  the  increase  in  gravel  content  with 
depth.  During  the  last  6000  years  sediment  deposition  has  likely  been  the  result 
of  material  being  transported  from  the  north  and  northeast.  As  there  appears  to 
be  few  modem  reef  organisms  producing  sediment  in  the  area  the  existing 
sediments  will  be  degrading  and  thus  will  become  fine  grained  through  time.  It  is 
to  be  expected  that  the  finer  material  will  be  deposited  in  hydrodynamically  quiet 
areas  such  as  small  basins  or  areas  protected  by  bathymetric  highs.  This  may 
account  for  sediments  becoming  increasingly  finer  grained  at  the  western  edge  of 
the  northern  reef  flat  which  is  protected  to  the  east  by  the  platform  situated 
between  the  central  and  northern  reef  flat.  The  coarse  material  at  the  north  end 
of  the  central  reef  flat  may  be  the  result  of  being  located  on  a  topographic  high. 
Also,  the  actual  boundary  of  the  reef  flat  as  determined  by  the  sidescan  sonar  may 
not  be  accurate.  The  relatively  high  percentage  of  gravel  in  the  central  portion 
of  the  survey  area  appears  to  correlate  with  the  central  axis  of  the  northeast 
trending  trough  and  may  reflect  higher  energy  conditions  due  to  tide-induced 
circulation  patterns. 

In  summary  the  surficial  geology  can  be  broadly  subdivided  into  reef  flats  and 
shallow  shelf  sediments.  The  shelf  sediments  appear  to  become  coarser  with 
subsurface  depth  especially  to  the  northwest.  The  shelf  sediments  show  relatively 
little  variability  except  for  coarser  grained  material  at  the  northern  end  of  the 
central  reef  flat,  the  axis  of  the  northeast  trending  trough,  and  to  the  west  of  the 
survey  area.  The  sampling  bias  prohibits  a  more  rigorous  classification. 
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2.5  Conclusions 

Data  from  a  total  of  58  sample  sites  augmented  with  sidescan  sonar  data  and  video 

imagery  provides  the  following  interpretation  of  the  surficial  geology  of  the 

survey  area: 

1 .  The  bathymetry  shows  a  broad  shallowing  trend  to  the  northwest.  Total 
relief  in  the  area  is  approximately  11m.  Bathymetric  features  include  a 
northeast  trending  trough  to  the  west  and  a  series  of  pinnacles,  3m  and  4m 
in  height,  in  the  southern  portion  of  the  survey  area. 

2.  Sediment  textural  analysis  indicates  that  the  majority  of  sediments  can  be 
classed  as  silts  and  clayey  sUts  except  for  samples  located  at  the  east  central 
portion  of  the  survey  area  and  along  the  axis  of  the  northeast  trending 
trough  to  the  west  of  the  survey  area. 

3.  There  appears  to  be  an  increase  in  grainsize  with  subsurface  depth, 
especially  in  the  northwest, 

4.  The  surficial  geology  can  be  classified  as  two  distinct  units:  reef  flat  and 
shallow  shelf  sediments.  While  there  appears  to  be  some  variability  in 
grainsize  the  patchy  nature  of  the  sampling  site  distribution  prohibits  a 
more  rigorous  classification. 
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3 .  Classification  Processing 

3.1  Procedure 

An  acoustic  pulse  impinging  on  a  seabed  is  reflected  and  scattered  by  both  the 
interface  between  the  water  and  the  seabed  and  from  volume  inclusions  within  the 
seabed.  The  returned  echo  pulse  is  shaped  by  the  seabed  roughness  and 
geoacoustic  properties.  We  make  use  of  the  echo  shape  in  order  to  perform  a 
seabed  classification. 

First,  we  generate  a  great  number  of  shape  parameters,  currently  166,  from  each 
echo.  We  cannot  easily  visualize  or  process  so  great  a  number  of  parameters,  but 
fortunately  most  parameters  carry  limited  information  or  redundant  information. 
This  means  that  for  a  collection  of  echoes  the  variance  of  any  individual 
parameter  is  small  (limited  information)  or  the  covariance  of  any  individual 
parameter  with  the  other  parameters  is  small  (redundant  information).  Second, 
we  use  principal  component  analysis  to  determine  which  parameters  carry 
information  and  extract  these  parameters.  Finally,  we  analyze  these  extracted 
parameters  using  either  supervised  or  unsupervised  classification  techniques. 

Preprocessing  is  performed  prior  to  echo  shape  analysis.  First,  the  seabed  echo  is 
located  (bottom  pick).  Next,  an  ensemble  of  consecutive  echoes  is  assembled. 

Each  echo  in  the  ensemble  is  aligned,  an  average  computed,  and  a  single  output 
ping  generated.  Finally,  echo  pulse  contractions  or  dilations  caused  by  altitude 
changes  are  removed  from  the  average  ping. 

Features  related  to  the  shape  of  the  echo  are  extracted  from  the  preprocessed 
ping.  We  call  this  collection  of  features  the  full  feature  vector.  Our  current 
implementation  produces  166  feature  vector  elements  from  five  algorithms:  a 
histogram  of  the  distribution  of  the  amplitudes  in  the  echo;  quantiles  of  the 
distribution  of  the  amplitudes  in  the  echo;  integrals  of  the  amplitudes  to  various 
times  in  the  echo  and  ratios  of  these  integrals;  Fourier  spectrum  amplitude 
coefficients;  and  wavelet  coefficients. 

Full  feature  vectors  from  a  variety  of  seabed  types  are  then  used  in  a  principal 
component  analysis.  A  covariance  matrix  of  dimension  166  x  166  is  produced  and 
the  eigenvalues  and  eigenvectors  are  calculated.  Eigenvectors  are  a  set  of 
orthonormal  basis  vectors  spanning  the  covariance  matrix  which  can  be  used  in 
conjunction  with  the  eigenvalues  to  account  for  the  energy  in  the  covariance 
matrix  most  rapidly  in  a  least- squares  sense.  We  have  determined  that  with  only  3 
out  of  166  eigenvectors  we  can  typically  account  for  over  95%  of  the  covariance 
produced  from  several  thousand  pings  spanning  a  wide  variety  of  seabed  types. 
We  therefore  use  the  eigenvectors  corresponding  to  the  three  largest  eigenvalues 
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as  weights  to  reduce  the  166  Ml  feature  vector  elements  to  a  reduced  feature 
vector  of  3  elements  -  features  1,  2,  and  3.  Features  1,  2,  and  3  are  the  elements 
of  the  reduced  feature  vectors  and  correspond  to  the  directions  associated  with 
the  largest,  second  largest  and  third  largest  eigenvalues  respectively. 

The  principal  component  analysis  provides  us  with  a  set  of  reduced  feature 
vectors  which  contain  most  of  the  covariance  energy.  Since  the  covariance  within 
the  same  seabed  type  is  less  than  the  covariance  between  different  seabed  types, 
we  expect  the  reduced  feature  vectors  to  be  clustered  around  locations  in  reduced 
feature  space  corresponding  to  a  seabed  type.  Two  general  methods  of 
determining  the  statistics  (mean  and  covariance)  of  the  locations  around  which 
clusters  occur  are  available.  Note  that  the  covariance  of  a  cluster  describes  the 
orientation  and  size  of  a  cluster  in  reduced  feature  space  and  is  not  related  to  the 
covariance  between  pings  used  by  the  principal  component  analysis.  First,  the 
statistics  of  each  set  of  reduced  feature  vectors  for  each  seabed  type  can  be 
calculated  (supervised  classification);  or  second,  an  assumption  about  the 
distribution  of  points  within  a  cluster  can  be  made  and  then  the  location,  size,  and 
number  of  clusters  can  be  adjusted  and  a  chi-squared  statistic  can  be  used  to 
determine  whether  the  adjustment  has  resulted  in  a  better  or  poorer  fit  to  the 
assumed  distribution  (unsupervised  classification).  Supervised  classification  is 
most  appropriate  if  ground  truthing  exists  and  known  “acoustic  signatures  can  be 
assigned  to  each  seabed  type.  Unsupervised  classification  is  most  appropriate 
where  there  is  little  or  no  ground  truthing  and  therefore  no  representative  echoes 
from  a  seabed  type  can  be  located  and  used  to  determine  the  acoustic  signature.  In 
this  case,  the  clustering  analysis  will  identify  acoustic  regimes  -  areas  which  have 
the  same  echo  shape  characteristics  -  thus  allowing  classification  into  areas  of  like 
and  unlike  echo  shapes. 

Once  the  statistics  of  the  clusters  corresponding  to  the  various  seabed  types  or 
acoustic  regimes  have  been  determined,  classification  is  made  by  determining  to 
which  cluster  any  new  ping’s  reduced  feature  vector  most  likely  belongs. 


3.2  Data  Acquisition 

An  ISAH-S  acoustic  data  acquisition  system  was  used  to  log  the  following  data: 

1 .  Simrad  EA300P  200  kHz  echo  sounder  waveforms 

2.  Raytheon  DSF6000  24  kHz  echo  sounder  waveforms 

The  waveform  data  were  obtained  from  a  post-detection  point  in  the  echo 
sounders  and  therefore  represent  the  envelope  of  the  time-varying  gain  (TVG) 
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corrected  signal  at  the  transducer.  The  transducer  specifications  are  not  known 
but  are  believed  to  be  in  the  range  of  10°  beamwidth  for  the  200  kHz  sounder  and 
20°  beamwidth  for  the  24  kHz  sounder.  Pulse  lengths  were  0.1ms  (200  ffiz)  and 
0.2ms  (24  kHz).  The  waveform  envelopes  were  sampled  at  20  kHz  providing  a 
possible  bandwidth  to  Nyquist  of  10  kHz.  Actual  data  bandwidth  is  only  several 
kHz  in  both  frequency  cases.  Data  were  logged  to  Exabyte  8mm  helical  scan  tape. 
Examples  of  200  kHz  and  24  kHz  waveforms  are  shown  in  Figures  3.1  and  3.2. 

An  ISAH-HYDAS  hydrographic  data  acquisition  system  was  used  to  log  the  GPS 
position  (GPS  manufacturer  is  unknown).  Unfortunately  the  differential  link 
operation  was  intermittent  so  the  positioning  is  only  generally  accurate  to  SA 
performance  (<100  m).  Data  were  logged  to  quarter  inch  cartridge  (QIC) 
DC6150  streamer  tapes. 
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Figure  3.1  200  kHz  waveforms. 
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3.3  Data  Pre-processing 

Positions  were  processed  using  QTC's  HYPS  processing  software.  This  involved 

downloading  the  raw  GGA  records  from  the  streamer  tapes,  editing  out 

obviously  bad  points,  and  filtering  the  remaining  positions  with  a  Kalman  filter. 

HYPS  outputed  a  position/time  table  in  ASCII  format. 

The  acoustic  waveforms  were  prepared  for  seabed  classification  by  applying  the 

following  processing  steps: 

1  The  waveforms  (traces)  were  downloaded  from  8mm  tape  to  disk  for  more 
rapid  access.  Over  the  survey  area  more  that  44,000  (200  kHz)  and  70,000 
(24  kHz)  waveforms  were  acquired. 

2.  The  seabed  (bottom)  was  tracked  and  the  water  depth  used  to  produce  the 
relative  bathymetry  figure. 

3.  A  5-fold  stack  (averaging)  of  the  traces  was  performed.  Prior  to 
averaging  the  bottom  picks  of  all  traces  are  aligned  to  reduce  smearing 
introduced  by  ping-to-ping  depth  variations.  The  ping  rates  were  1.1  Hz 
(200  kHz)  and  1.8  Hz  (24  kHz)  and  vessel  speed  was  approximately  2  m/s. 
Therefore  the  5-fold  stack  represents  a  horizontal  separation  of 
approximately  9m  at  200  kHz  and  5.6m  at  24  kHz.  The  stacking  was  not 
rolled  along  (moving  average)  therefore  the  data  were  spatially  decimated: 
every  5-fold  stack  was  output  at  200  kHz;  every  other  5-fold  stack  was 
output  at  24  kHz.  Total  number  of  traces  were  therefore  8800  at  200  kHz 
and  7000  at  24  kHz.  These  trace  densities  represent  one  averaged  trace  at 
approximately  every  10m  along  track.  Water  depth  below  transducer  for 
the  survey  area  averages  about  25m;  this  means  that  the  acoustic  footprint 
sizes  (-3dB)  are  about  4m  at  200  kHz  and  9m  at  24  kHz.  The  ou^ut  trace 
spacing  of  1  trace  per  10m  represents  only  1-2  times  the  acoustic  footprint 
size  and  therefore  does  not  seriously  degrade  the  possible  spatial  resolution. 

4  Depth  effects  are  removed  by  normalizing  all  echoes  to  a  constant  virtual 
depth.  The  differential  depth  effect  is  not  large  for  these  data  as  water 
depth  is  reasonably  constant. 

5_  The  echo  shape  is  analyzed  using  the  5  algorithms  described  in  Section  3.1 
above.  Each  echo  is  characterized  by  166  features. 
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6.  Principal  component  analysis  is  used  to  reduce  the  166  full  feature  vector 
to  3  principal  components,  the  reduced  feature  vector  consisting  of  Ql,  Q2, 
and  Q3. 


3.4  Unsupervised  Classification 

Unsupervised  classification  means  that  the  classes  are  automatically  chosen  by  the 
system  using  a  Chi-squared  criteria.  For  the  unsupervised  classification  the 
principal  component  analysis  was  run  on  every  10th  trace  processed  by  steps  1-5 
in  Section  3.3  above.  The  24  kHz  and  200  kHz  reduced  features  were  analyzed 
separately. 

The  unsupervised  clustering  program  assumes  that  the  reduced  feature  vectors 
cluster  around  the  various  bottom  types.  Each  bottom  type  therefore  can  be 
represented  by  a  location  in  reduced  feature  vector  space  (Ql,  Q2,  Q3  space)  and 
a  hyperellipse  representing  the  1  standard  deviation  surface  of  the  cluster.  If  we 
assume  that  the  distribution  of  points  within  a  cluster  is  Gaussian,  we  can  adjust 
the  size  and  orientation  of  the  hyperellipse  enclosing  the  points  until  the  Chi- 
squared  fit  to  a  Gaussian  is  minimized.  If  the  Chi-squared  fit  cannot  be  reduced 
to  the  number  of  degrees  of  freedom  minus  1,  in  this  case  1,  the  cluster  is  split 
into  two  clusters  and  the  process  continued. 

In  this  way  both  the  200  kHz  and  24  kHz  can  be  clustered  into  4  different  classes. 
By  inspection  of  the  location  of  these  4  classes  with  the  ground  truing  described 
in  Section  2  above,  the  following  class  descriptions  were  chosen:  silt,  sand/gravel, 
reef  1,  and  reef  2.  Reduced  feature  vectors  in  Ql,  Q2,  Q3  space  and  associated 
classes  for  the  two  different  frequencies  are  plotted  in  Figure  3.3  and  3.4. 

3.5  Supervised  Classification 

Supervised  classification  defines  classes  based  on  areas  of  a  known  bottom  type. 
Feature  vectors  from  these  known  areas  are  generated,  reduced  using  principal 
component  analysis,  and  the  class  descriptions  are  calculated  for  each  known 
bottom  type.  Unfortunately,  ground  truthing  for  the  survey  area  was  insufficient 
to  perform  a  good  supervised  classification.  The  original  intent  was  to  augment 
the  unsupervised  classification  classes  with  any  additional  known  classes; 
however,  the  ground  truthing  quality  on  the  survey  lines  is  inadequate  to  do  t\As. 
The  positioning  of  the  acoustic  and  ground  truthing  data  must  be  of  high  quality 
and  must  coincide  and  the  ground  truthing  itself  must  be  accurate  to  allow 
supervised  classification.  As  processing  of  the  ground  truthing  from  the  survey 
area  becomes  available  it  may  be  possible  to  augment  the  unsupervised  results. 
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Figure  3.3  200  kHz  Cluster  Plot. 
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Figure  3.4  24  kHz  Cluster  Plot. 
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4 .  Classifications  Results 

Maps  showing  the  classification  results  are  presented  in  this  section.  Discussion 
of  results  is  deferred  to  Section  5.  Two  types  of  presentations  are  made  for  each 
frequency:  first,  trackplot  maps  are  shown  with  track  location  colour  coded  by 
bottom  type;  and  second,  proportional  composition  plots  for  each  frequency  are 
presented.  This  section  presents  proportional  composition  plots  of  the  entire 
survey  area  on  one  plot;  proportional  composition  plots  made  on  a  trackline  basis 
are  shown  in  Appendix  B. 

Two  types  of  trackplots  are  shown  for  each  frequency.  The  first  trackplot  shows 
all  four  bottom  types  as  different  colours,  the  second  trackplot  plots  the 
sand/gravel,  reef  1,  and  reef  2  classes  as  all  the  same  colour.  The  two  colour 
trackplots  show  the  effect  of  echo  sounder  frequency  more  clearly  than  the  four 
colour  plots. 

The  proportional  composition  plots  are  generated  by  taking  an  ensemble  of  20 
successive  classifications  and  outputting  the  number  of  each  class  within  the 
ensemble.  These  plots  are  useful  because  they  show  that  the  seabed  is  not 
generally  homogenous  and  show  how  transitions  between  various  bottom  types 
occur.  The  individual  trackline  plots  shown  in  Appendix  B  are  generated  from 
ensembles  of  10  traces. 

4.1  Classifications 

Figures  4.1  -  4.4  show  the  trackplots  described  above.  Figures  4.5  -  4.6  show  the 
proportional  compositions  of  the  two  frequencies. 

These  results  are  discussed  in  Section  5. 
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Figure  4. 1  Classification  of  200  kHz  echo  sounder  data. 
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Figure  4.2  Classification  of  24  kHz  echo  sounder  data. 
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Figure  4.3  Classification  of  200  kHz  echo  sounder  data.  Classes  sand/gravel, 
reef  1,  and  reef  2  are  all  plotted  as  yellow. 
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Figure  4.4  Classification  of  24  kHz  echo  sounder  data.  Classes  sand/gravel,  reef 
1,  and  reef  2  are  all  plotted  as  yellow. 
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Sand/Gravel  Proportional  Composition  Plot  of  200  kHz 
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5.  Discussion  and  Conclusions 

5.1  General  Results 

The  cluster  analysis  produced  four  different  classes  in  both  the  24  kHz  and  200 
kHz  cases.  The  clusters  and  results  of  the  unsupervised  classification  are  shown 
in  Figures  3.3  and  3.4.  The  general  orientation  of  classes  is  the  same  for  both 
frequencies  although  the  200  kHz  data  exhibit  a  much  greater  spread  of  the  4th 
class  (plotted  as  cyan)  than  does  the  24  kHz  data.  Inspection  of  the  locations 
associated  with  samples  from  each  cluster  class  leads  to  classification  into  the  four 
categories:  silt,  sand/gravel,  reef  1  and  reef  2.  The  classes  and  their  orientation 
in  reduced  feature  vector  space  indicate  a  coarsening  of  sediment.  Further 
ground  truthing  is  required  to  refine  these  classes  and  confirm  their  naming. 

Actual  classification  results  are  presented  in  Figures  4.1  and  4.2.  Each  plotted 
classification  represents  the  bottom  type  over  about  10m;  the  size  of  the  colour 
dots  in  the  figures  also  represent  about  10m.  There  is  very  good  general 
agreement  between  the  results  fi:om  the  two  frequencies,  although  the  24  kHz 
result  shows  more  coarse  material.  This  effect  is  discussed  in  Section  5.2  below. 

The  proportional  composition  plots  (Figures  4.5  and  4.6)  show  how  the  seabed 
grades  from  one  bottom  type  to  another.  Figures  4.5  and  4.6  were  made  by 
collecting  the  classification  results  for  all  survey  lines  including  turns  and  plotting 
them  sequentially  (basically  replaying  the  survey  in  proper  time  sequence). 
Proportional  composition  plots  made  on  a  line-by-line  basis  (waypoint  to 
waypoint)  are  presented  in  Appendix  B.  These  plots  are  most  effectively 
interpreted  by  comparing  the  24  kHz  to  the  200  kHz  plot.  In  the  200  kHz  plot  we 
generally  see  a  greater  proportion  of  silt,  and  that  some  areas  are  exclusively  silt. 
In  the  24  kHz  plot  much  of  the  silt  class  has  moved  into  the  sand/gravel  class. 

The  reef  1  class  is  also  larger  on  the  24  kHz  plot;  the  reef  2  class  seems  to  occur 
about  as  frequently  as  in  the  200  kHz  case  but  the  proportion  is  slightly  less. 

These  effects  are  discussed  in  Section  5.2  below. 


5.2  Effect  of  Echo  Sounder  Frequency 

The  effect  of  the  echo  sounder  frequency  is  most  clearly  seen  in  Figures  4.3  and 
4.4.  In  these  both  these  figures  the  silt  class  is  plotted  in  red  and  all  other  classes 
(not  silt)  plotted  in  yellow.  We  see  on  the  200  kHz  plot  (Figure  4.3)  the  only 
substantid  areas  of  yellow  are  the  three  reef  flat  areas  described  in  Section  2. 

This  result  is  consistent  with  the  relatively  clean  reef  flat  areas  shown  in  the  video 
and  interpreted  from  the  general  geology  of  the  area.  The  24  kHz  plot  shows 
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much  more  yellow  (non-silt)  areas.  The  only  significant  silt  areas  bracket  the 
central  reef  flat;  the  areas  over  the  northern  and  southern  reef  flats  show  much 
less  silt  that  on  the  200  kHz  plot.  There  is  no  silt  area  on  the  24  kHz  plot  which  is 
not  silt  on  the  200  kHz  plot  except  for  several  isolated  points. 

Interestingly,  the  24  kHz  plot  shows  the  non-silt  class  gradually  increases  towards 
the  northwest.  Table  2.2  and  Figure  2.1  show  the  locations  for  which  gravel 
content  with  depth  cores  exist  and  the  increase  in  gravel  content  with  depth.  The 
cores  generally  show  an  increase  in  gravel  content  with  depth,  i.e.,  the  bottom 
consists  of  a  gravel  bottom  overlain  by  silt.  It  is  likely  that  at  200  kHz  there  is  no 
penetration  and  hence  the  seabed  is  classified  as  silt;  whereas  at  24  kHz  the 
increased  penetration  into  the  bottom  allows  the  underlying  sediment  to  be 
classified. 

The  proportional  composition  plots  (Figures  4.5  and  4.6)  show  much  the  same 
result.  The  24  kHz  echo  sounder  seems  to  be  able  to  penetrate  much  of  the  silt 
overlying  the  sand/gravel  and  reef  1  classes  and  classify  the  underlying  substrate. 
The  reef  2  class,  perhaps  associated  with  newer  reef  material  and  associated  with 
the  pinnacles  does  not  show  an  increase  in  proportion.  This  indicates  that  these 
pinnacles  are  classified  as  the  same  in  both  the  200  kHz  and  24  kHz  cases, 
probably  because  the  pinnacles  do  not  contain  pockets  of  silt  as  do  the  reef  flats. 

Unfortunately  the  core  locations  are  generally  clustered  in  the  northwest  quadrant 
and  are  not  on  the  lines  surveyed  with  the  echo  sounders.  Additional  cores  on  the 
survey  lines  and  over  the  various  bottom  types  shown  on  the  two  plots  (silt  on 
200  kHz  and  24  kHz;  not  silt  on  200  kHz  and  24  kHz;  sdt  on  200  kHz  and  not  silt 
on  24  kHz)  would  be  quite  useful. 


5.3  Accuracy 

In  section  2.5  conclusions  regarding  the  ground  truthing  interpretation  were 
made.  These  points  are  reiterated  and  examined  in  light  of  the  acoustic 
classification: 


1.  The  bathymetry  shows  a  broad  shallowing  trend  to  the  northwest.  Total 
relief  in  the  area  is  approximately  11m.  Bathymetric  features  include  a 
northeast  trending  trough  to  the  west  and  a  series  of  pinnacles,  3m  and  4m 
in  height,  in  the  southern  portion  of  the  survey  area. 


The  pinnacles  are  classified  as  reef  2  on  both  the  24  kHz  and  200  kHz  plots 
(Figures  4.1  and  4.2).  The  increase  in  sediment  coarseness,  especially  at  24 
kHz  is  also  consistent  with  the  shallowing  trend  as  fine  sediments  such  as 
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silt  tend  to  be  depositied  in  deeper  basins  and  other  areas  protected  from 
currents  and  wave  action. 

2.  Sediment  textural  analysis  indicates  that  the  majority  of  sediments  can  be 
classed  as  silts  and  clayey  silts  except  for  samples  located  at  the  east  central 
portion  of  the  survey  area  and  along  the  axis  of  the  northeast  trending 
trough  to  the  west  of  the  survey  area. 

Both  the  24  kHz  and  200  kHz  datasets  could  be  divided  into  only  4  different 
classes.  Based  on  the  above  point  we  have  named  these  classes  silt, 
sand/gravel,  reef  1  and  reef  2.  We  suspect  reef  1  is  old  reef  and  reef  2  is 
newer  reef  associated  with  the  several  piimacles  found  on  the  old  reef  flats. 

3.  There  appears  to  be  an  increase  in  grainsize  with  subsurface  depth, 
especially  in  the  northwest. 

The  increase  in  the  non-silt  class  towards  the  northwest  for  the  24  kHz  data 
relative  to  the  200  kHz  data  indicates  that  as  the  silt  thickness  decreases  the 
24  kHz  is  able  to  penetrate  the  silt  and  classify  the  underlying  substrate 
whereas  the  200  kHz  never  penetrates  into  the  substrate  and  therefore 
always  classifies  the  overlying  silt. 

4.  The  surficial  geology  can  be  classified  as  two  distinct  units:  reef  flat  and 
shallow  shelf  sediments.  While  there  appears  to  be  some  variability  in 
grainsize  the  patchy  nature  of  the  sampling  site  distribution  prohibits  a 
more  rigorous  classification. 

The  reef  flat/non  reef  flat  areas  have  been  identified  in  both  the  200  kHz 
and  24  kHz  data. 

In  total,  agreement  with  known  ground  truthing  interpretation  is  very  good. 

Individual,  detailed  comparisons  with  the  actual  ground  truthing  data  is 

impossible  as  the  positioning  of  the  acoustic  data  is  not  good  enough  and  most 

sample  sites  are  not  on  the  lines  of  the  acoustic  survey. 


5.4  Conclusions 

The  purpose  of  this  report  is  to  present  results  of  seabed  classification  processing 
using  the  QTC  CVIEW  system  on  the  Florida  Keys  data  and  to  compare  with 
known  ground  truthing.  We  feel  that  the  classification  is  good  for  several 
reasons: 
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1 .  There  is  very  good  agreement  between  the  200  kHz  and  24  kHz  datasets, 
especially  over  the  reef  2  areas. 

2.  Areas  of  differences  between  classifications  at  the  two  frequencies  show  an 
increased  sediment  coarseness  in  the  24  kHz  results.  This  is  consistent  with 
the  greater  penetration  expected  at  the  lower  frequency  and  the  cores 
which  show  an  increasing  sediment  coarseness  with  depth. 

3.  Generally,  the  QTC  VIEW  classification  results  are  consistent  with  the 
known  geology  and  ground  truthing  of  the  area.  We  know  that  there  is 
little  variation  in  seabed  type  from  the  ground  truthing;  this  lack  of 
variation  (only  4  classes  and  not  much  along  track  variability  in  the 
acoustic  results)  is  exhibited  by  the  acoustic  data. 

4.  The  4  classes  were  generated  totally  automatically  by  using  a  Chi-squared 
clustering  program.  Both  the  200  kHz  and  24  kHz  data  produced  4  classes; 
for  both  frequencies  the  classes  were  ordered  the  same  way.  These  facts 
indicate  that  the  echo  sounders  acoustically  see  the  seabed  in  a  similar 
fashion  and  give  us  some  confidence  that  these  results  can  be  extended  to 
other  areas.  Supervised  classification  with  a-priori  knowledge  can  be  very 
effective,  especially  within  the  area  of  the  supervised  inputs,  but  can  prove 
to  be  much  more  difficult  to  extend  to  other  areas. 


A  major  area  for  additional  seabed  classification  processing  would  be  to  use  a 
multi-spectral  approach  to  seabed  classification.  In  this  approach  the  24  kHz  and 
200  kHz  feature  vectors  are  combined  and  principal  component  and  cluster 
analysis  performed  on  the  combined  feature  vectors.  This  should  yield  a 
combined  seabed/upper  sub-bottom  classification.  We  believe  that  this  approach 
can  be  extended  to  the  deeper  sub-bottom  by  adding  an  even  lower  frequency  to 
the  24  kHz  and  200  kHz  data.  Additional  studies  measuring  the  actual  penetration 
into  the  sub-bottom  at  the  various  frequencies  would  aid  the  interpretation  greatly 
as  the  hypothesis  that  the  24  kHz  is  classifying  the  substrate  could  be  confirmed. 

The  acoustic  classification  can  be  interpreted  only  in  a  general  sense  because  the 
ground  truthing  and  acoustic  survey  positioning  is  of  only  a  general  nature.  A 
more  rigorous  program  combining  detailed  positioning  and  ground  truthing  with 
an  accurate  acoustic  survey  is  a  necessary  next  step. 

If  this  experiment  were  to  be  repeated  or  extended,  we  would  make  several 
recommendations; 


1 .  The  ground  truthing  program  should  be  more  closely  coordinated  with  the 
acoustic  program.  This  would  allow  interpretations  and  classification 
success  estimates  to  be  made  on  a  more  local  scale  and  with  greater 
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confidence.  This  would  also  allow  supervised  classification  results  to  be 
included  with  the  unsupervised  classification  results. 

2.  Surficial  grab  samples  and  visual  seabed  identifications  should  be  made 
frequently  over  as  wide  a  variety  of  bottom  types  as  possible.  This  would 
allow  the  limits  of  classification  resolution  to  be  more  properly 
understood.  This  would  also  allow  supervised  classification  to  augment  the 
results. 

3.  A  multi-frequency  experiment  combining  several  echo  sounder  frequencies 
and  an  experiment  to  determine  depth  of  penetration  at  the  various 
frequencies  should  be  performed.  This  experiment  could  be  used  as  to 
investigate  the  possibility  of  extending  the  seabed  classification  to  sub¬ 
bottom  classification. 
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6 .  Appendix  A:  Ground  Truthing  Data 

Legend  to  Ground  Truthing  Data  supplied  by  EDRD; 

Core  logs  from  vessel  Planet:  Identified  by  name  Planet  in  cruise  location.  Site 
number  used  in  the  Figure  2.2  and  Table  3.1  is  the  numeric  part  of  the  station 
number. 

Core  logs  from  vessel  Pelican:  Identified  by  table  title  kw-pe-gc-###.  Site 
number  used  in  Figure  2.2  and  Table  3.1  is  the  numeric  part  of  the  title. 
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7.  Appendix  B:  Proportional  Composition  Plots 

Figure  B.l  Waypoint  locations. 
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Figure  B.4  WP  48  -  WP  47. 
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Figure  B.5  WP46-WP45. 
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Figure  B.7  WP42-WP41. 


Sand/Gravel  Proportional  Composition  Plot  of  200  kHz  Run007 
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Figure  B.8  WP  68  -  WP  69. 
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Figure  B.9  WP  66  -  WP  67. 

Sand/Gravel  Proportional  Composition  Plot  of  200  kHz  Run009 
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Figure  B.ll  WP  62  -  WP  63. 


Sand/Gravel  Proportional  Composition  Plot  of  200  kHz  Run028 
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